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Abstract. Irregular arrangements of vesicles rilled with excitable and pre- 
cipitating chemical systems are imitated by Voronoi automata — finite-state 
J.^ machines defined on a planar Voronoi diagram. Every Voronoi cell takes four 

states: resting, excited, refractory and precipitate. A resting cell excites if it 
has at least one excited neighbour; the cell precipitates if a ratio of excited 
cells in its neighbourhood to its number of neighbours exceed certain thresh- 
old. To approximate a Voronoi diagram on Voronoi automata we project a 
planar set onto automaton lattice, thus cells corresponding to data-points are 
excited. Excitation waves propagate across the Voronoi automaton, interact 
with each other and form precipitate in result of the interaction. Configuration 
of precipitate represents edges of approximated Voronoi diagram. We discover 
^ relation between quality of Voronoi diagram approximation and precipitation 

q threshold, and demonstrate feasibility of our model in approximation Voronoi 

I I diagram of arbitrary-shaped objects and a skeleton of a planar shape. 

' ^ Keywords: Voronoi diagram, automata, discrete networks, wave dynamics 

o 

o 

1. Introduction 

A regular, or irregular but manually designed, arrangements of vesicles filled 
with excitable chemical mixtures bear huge computational potential [B]. When 
vesicles are in close, at least in a diffusion terms, contact with each other, via tiny 
pores, excitation waves can pass from one vesicle to its close neighbour. Excitation 
wave-fragments keep their shape, more or less constant, inside each BZ-vesicle. A 
; i wave-fragment passing from one vesicle to another it contracts, due to the restricted 

size of the connecting pore. When two or more wave-fragments collide inside a vesi- 
cle they can annihilate, deviate, or multiply. When interpreting presence/absence 
of wave-fragments in any given as a value of Boolean variable we can implement 
all basic operations of a Boolean logic via collisions between wave-fragments in a 
vesicle. In computer experiments we designed a binary adder in a hexagonal ar- 
ray of vesicles filled with excitable chemical mixture [B] , built polymorphic logical 
gates (switching between XNOR and nor) by using illumination to control outcomes 
of inter-fragment collisions [7] and geometry-modulated complex arithmetical cir- 
cuits |27j . Our theoretical ideas and results got experimental chemical laboratory 
back up - results on information transfer between Belousov-Zhabotinsky mixture 
enclosed in lipid membrane are successful [231 I33j • 
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(c) 

Figure 1 . Examples of Voronoi representation of bubble conglom- 
erates: (a) rhodamin, vegetable oil, water and silicon oil (poly- 
dimethylsiloxane), (b) blue food coloring and vegetable oil, (c) rho- 
damin, water and silicon oil mixtures. 



While in computer models regular arrangement of uniform vesicles is effort- 
less real-life experiments bring nasty surprises. Usually vesicles are different sizes, 
they do not form a hexagonal lattice as a rule, and they may be unstable, coa- 
lescence transforms fine-grained networks of elementary vesicle-processors into a 
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coarse-grained assembles of monstrous vesicular structures. What kind of com- 
putation can be done on an irregular arrangement of non-uniform vesicles? We 
address the question by representing the vesicle assembles by automata networks 
and studying how planar subdivision, Voronoi diagram, problem can be solved in 
such vesicle-automata. 

We abstract vesicle assembles as planar Voronoi diagrams (Fig.[lJ of planar sets, 
points of which arc centres of the vesicles. Voronoi diagram is routinely as approx- 
imation of arrangements of discs [55] and sphere packing [3TJ [5p , 32 . The diagram 
is also used in structural analysis of liquids and gases [9] , and protein structure |36j , 
and to model dense gels [JT] and inter-atomic bonds [2E]- The Voronoi diagrams 
are introduced in Sect. [2] We assume that every cell of a Voronoi diagram is a 
finite-state that takes four states and updates its states depending on states of its 
first and second order neighrbours. We design a cell-state transition function which 
combines generalised, and highly-abstracted, properties of both excitable and pre- 
cipitating chemical media: a local disturbance gigves birth to quasi-circular waves 
of excitation while collisions between the waves lead to precipitation. The Voronoi 
automaton is defined in Sect. [3j 

The problem solved by Voronoi automats is the approximation of Voronoi dia- 
gram. Approximated Voronoi diagram is much more coarse-grained than Voronoi 
diagram on which excitable- precipiating automaton is built. To approximate a 
Voronoi diagram on Voronoi automata we project a planar set onto automaton 
lattice, thus cells corresponding to data-points are excited. Excitation waves prop- 
agate across the Voronoi automaton, interact with each other and form precipitate 
in result of the interaction. Configuration of precipitate represents edges of approx- 
imated Voronoi diagram. In our model precipitation depends on a local density of 
excitation — precipitation threshold. For low precipitation threshold the medium 
becomes cluttered with meaningless clusters of precipitate, for high threshold few 
domain of precipitation is formed. Our quest for optimal threshold of precipiation 
is narrated in Sect. [4j In Sect. [5] we show optimal threshold found can be equally 
used to approximate Voronoi diagram of arbitrary planar shapes and also skeleton 
of planar shape. 

2. Voronoi diagram 

Let P be a nonempty finite set of planar points. A planar Voronoi diagram [10] 
of the set P is a partition of the plane into such regions, that for any element of 
P, a region corresponding to a unique point p contains all those points of the plane 
which are closer to p than to any other node of P. A unique region vor(p) = {z£ 
R 2 : d(p, z) < d(p, m)Vm 6 R 2 , m ^ z} assigned to point p is called a Voronoi cell 
of the point p [37] ■ The boundary dvor(p) of the Voronoi cell of a point p is built of 
segments of bisectors separating pairs of geographically closest points of the given 
planar set P. A union of VD(P) = U p& pdvor(p) all boundaries of the Voronoi cells 
determines the planar Voronoi diagram |37| . A variety of Voronoi diagrams and 
algorithms of their construction can be found in [5SJ |M] . 

Approximation of Voronoi diagrams with propagating patterns is based on time- 
to-distance transformation: to approximation a bisector separating planar points p 
and q we initiate growing patterns at p and q. The pattern travel the same distance 
from the sites of origination before they meet each other, The loci where the waves 
meet indicate sites of the computed bisector [TJ [5] . Precipitating reaction-diffusion 
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chemical media are proved to be an ideal computing substrate for approximation 
of the planar Voronoi diagram [3§1 HS1 13] ■ A Voronoi diagram can be approximated 
in a two-reagent medium. One reagent a is saturated in the substrate, drops of 
another reagent j3 are applied to the sites corresponding to planar points to be 
separated by bisectors. The reagent (3 diffuses in the substrate and reacts with 
reagent a. Colored precipitate is produced in the reaction between a and ft. When 
two or more waves of diffusing a meet, no precipitate is formed [18] . Thus uncolored 
loci of the reaction-diffusion medium represent bisectors of the computed diagram. 
A range of chemical precipitating processors is designed and working prototypes 
are tested in laboratory conditions [3"9l IT51 IT71 ITS1 IP9] . 

3. Voronoi automata 

A Voronoi automaton is a tuple V = (V(P), Q, N, u, /), where P is a finite planar 
set, V(P) = {V(p) : p <E P}, Q is finite set, N is a set of natural numbers and 
u : V(P) V(P) fe is second-order neighbourhood, < k < |P|, and / : Q fc -> Q 
is a cell-state transition function. 

Excitable-precipitating Voronoi automata studied in present paper are specified 
as follows. Cell state set has for elements, Q = {o, +, — , Thus we assign three 
excitation-related states — resting (o), excited (+) and refractory (— ) — to cells, 
and one precipitate state #. Cells update their state in discrete time. A state of 
cell V(x) are time step t e N is denoted as V(x) t . All cells update their states in 
parallel using the same cell-state transition function. 

Let w(p), p G P, be a first-order neighbourhood of a Voronoi cell Vor(p) from 
V(P): w(p) = {V(q) € V(P) : dV(p) H dV(q)}, i.e. a set of Voronoi cells which 
have common edges with V(x). A second-order neighbourhood u(p) is a set of 
neighbours of first-order neighbours of V(p): u{p) — C\v(q)£w(v(p)) w (V{<l)- Exam- 
ples of Voronoi cell neighbourhoods are shown in Fig. [2^,. Distributions of numbers 
of immediate and second neighbours and sizes of second-order neighbourhood u are 
given in Fig[2jD. Predominant number, over 40%, of cells have six first-order (imme- 
diate) neighbours. Almost half of the cells have 13 or 14 second-order neighbours. 
For any cell size of second-order neighbourhood u is a sum of first-order neighbours 
and second-order neighbours. For Voronoi automata with high-density packing of 
15K cells, half of the cells have either 19, or 20 or 21 neigbours. If we consider 
only immediate neighbourhood then V can be, in principle, seen as a slightly dis- 
torted, hexagonal lattice. However, by adopting second-order neighbourhood we 
are moving to a random structure with higher than hexagonal lattice coordinate 
number (a node in hexagonal lattice has 18-node second-order neighbourhood while 
dominating neighbourhood sizes in V are 19, 20 and 21). 

Transition form excited to refractory state is unconditional, i.e. takes place with 
regards to states of a cell's neighbours. A resting cell excites if it has at least 
one excited neighbour. In this particular model we take a refractory state and 
precipitate states are absorbing: one a cell takes either of these two state it does 
not update its state any longer. 

Neighbourhood sizes may differ between Voronoi cells therefore we use a state 
transition function, where a cell updates its state depending on a relative excitaton 
in its neghbourhood. We assume that precipitation occurs in a resting cell when a 
ratio of excited neighbours to a total number of neighbours exceeds some threshold 
r) e [0, 1]. Let <j(V(x) r ) be a number of excited cells in the cell V(x)'s second-order 
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Number of neighbours of Voronoi cell 



(b) 



Figure 2. Structure and properties of neighbourhoods, (a) Ex- 
ample of Voronoi tesselation with second-order neighbourhood 
structures highlighted for four cells. In each case central cell p 
is filled with blue (black in gray-scale reproduction) colour, first- 
order neighbours of w(p) are filled with red (dark-gray) colour, 
and neighbours of w(p) are filled with green (light-gray) colour. 
For each blue (black) Voronoi cell a second-order neighbourhood 
is a set of red (dark-gray) and green (light gray) Voronoi cells, 
(b) Distribution of first order, or immediate, neighbours (dotted 
line), second order neighbours (dashed line) and size of second or- 
der neighbourhood (solid line) in Voronoi automata of 15K cells. 
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neighbourhood, a{V{xf) = Y.v(y)eu(v(x)) \{ v iv) '■ v ivf = +}l then a cell updates 
its state by the following rule: 

{#, if V(x) t = o and cf t (x)/v(x) > 77 
+, if V{x)* = o and a t (x)/u(x) > 1 
, i£V(x)* = + 
o, otherwise 

In computational experiments we construct set P by fill a disc-container of radius 
480 units with up to 15K points. The points are packed at random but there 
is at least 5 units distance between any two points. Voronoi diagram V — on 
which Voronoi automaton is constructed - is calculated by a classical sweepline 
algorithm [2"T] . 



4. Constructing Voronoi diagram on Voronoi automata 

Let B be a set of planar points on which a Voronoi diagram 03 on V, 77 = 0.4, is 
approximated. We project B onto V and excite cells of V which are closer than 9 
units to points of B. Excitation waves spread on V. The excitation waves collide 
and precipitation occurs nearby sites of the waves' collisions. Configuration of cells 
in precipitate state represents edges of 03(B. An example of excitation-precipitation 
dynamics in V which approximates 23(B), where B is a planar set of seven points, 
is shown in Fig. [3j A scheme of the Voronoi diagram computed is shown in Fig. [7^,. 

Even in a single example shown in Fig. [3] it is clear that V does not calculate 
93(B) precisely, it rather approximates it with discrete domains of precipitation, 
and also introduces some noise (sites of precipitation not coinciding with bisectors 
of Q3(V)). Only threshold of precipitation 77 can be varied in our model. How are a 
degree of approximation and noise to signal ratio depends on value of 77? To answer 
we undertook a series of computational experiments, illustrated - for some values 
of 77 — in Fig. [4] 

For a set B of seven points we constructed 'classical' planar Voronoi diagram 
03(B) (Fig.Qi). Then for 77 = 0.2 . . . 0.525 we excited V with B, waited till V reaches 
its stationary configuration, where all cells are in either precipitate or refractory 
state, and compared the configurations of V and image of 03(B). For each 77 we 
recorded a fraction p of precipitate-cells in a final configuration of V(rj) to a total 
number of precipitate-cells in the same configuration (Fig. [5^). 

We found that while density v, calculated as a ratio of a number of precipitate- 
cells in stationary configuration V(rj) to a number of precipitate-cells in configu- 
ration of V(0.2), exponentially decreases with increase of precipitation threshold 
77 (Fig. [5]}) degree of approximation p shows more intriguing behaviour (Fig. [5^,). 
From 77 = 0.2 till 77 = 0.41 'signal-to-noise' ratio p(rf) grows almost exponentially. 
Trend of ,0(77) undergoes S-type transition between 77 = 0.41 and 77 = 0.453 with one 
pick /j(0.44) and two cavities p(0.43) and p(0.452). As illustrated in Fig. [| V with 
precipitation threshold in the zone around 77 = 0.4, onset of 'strange' behaviour of 
p(r])), demonstrates best approximation of 03(B), as recognised by unaided eye. 

5. Arbitrary-shaped planar objects and contours 

The Voronoi automaton V copes well with data set O of planar finite shapes. 
To approximate diagram 03 of arbitrary planar shapes we project the objects of O 
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(c) t = 20 (d) t = 44 



Figure 3. Approximation of Voronoi diagram 93(B) on Voronoi 
automaton V, r\ = 0.4. Points of B are projected onto V at time 
step t = 0. Seven Voronoi cells of V are excited and generated 
quasi-circular excitation waves. Configuration of V at time step 
t = 1, when wave have just started to develop, is shown in (a). 
The waves propagate outwards their initial stimulation sites and 
covert Voronoi cells they are occupying into refractory states (be) . 
When two or more waves collide precipitation occurs. By 44th 
step of the automaton development excitation extincts but domains 
of precipitate-cells represent edges of the approximated Voronoi 
diagram V. Resting cells are blank, excited cells red (dark gray), 
refractory cells gray, and precipitate cells black. 
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(a) VD 



(b) ■q = 0.2 



(c) r t = 0.225 



(d) jj = 0.25 



T] = 0.275 



(f) r, = 0.3 



(g) r) = 0.325 



(h) n = 0.35 



r] = 0.375 



(j) V = 0.4 



(k) t] = 0.425 



(1) r) = 0.45 



(m) ri = 0.4755 



(n) r\ = 0.5 



(o) ■q = 0.525 



Figure 4. Quality of approximation depends on rj. (a) Voronoi di- 
agram 9J(B) of seven planar points set B constructed by Fortuna's 
sweepline algorithm, (b) (o) stationary configurations of precipi- 
tate in automaton V excited by B for different r/ = 0.2 .. . 0.525. 
In snapshots (b)-(o) black pixels symbolise Voronoi cells of V in 
precipitate state, cells in refractory states are blank. 
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(a) (b) 

Figure 5. Characterisation of approximation of 23(B) (Fig. [4^,) 
by automaton V: (a) p vs 77 and (b) v vs 77. Markers (rhombs and 
circles) show actual values of 1/(77) an d p(v)> while solid lines are 
moving average trend line purely for eye guidance. 



(Fig. |6p,) onto V in such a manner that at time t = all cells of V covered (coinciding 
with pixels of) by shapes from O becomes excited (Fig. ^p) . Waves of excitation 
travel across V and provoke precipitation during their collisions (Fig. |6j;-e) . When 
all cells take refractory or precipitate state the automaton's configuration becomes 
stationary. The spatial distribution of precipitate represents Voronoi diagram of 
data shapes (see Fig. [6^ and Fig.^). A precipitate-enhancement of original position 
of date shape is a byproduct of Vs development. 

A skeleton of a planar contour is a set of centers of bitangent circles lying inside 
the contour |15) . Blum's grass-fire algorithm for computing skeleton employs prop- 
agating patterns (T3J 03J [TS] : to compute a skeleton we set a contour on fire and 
let the fire spread, wuench points where the advancing fire-fronts collide represent 
the skeleton. Most known algorithms of skeletonisation are based on Blum's ap- 
proach: simulation of grass- fire [30] , distance transform |38j , analytical construction 
of medial axis and topological thinning [551 HD] ■ 

To approximate a skeleton of a shape S (Fig. [8^) we project S onto V. Voronoi 
cells of V corresponding to black pixels of S are excited (Fig. [8Jd). Excitation 
waves propagate inside and outside of the domain S of initial excitation; we can 
ignore behaviour of outward waves. Waves travelling inside the contour S trigger 
precipitation at the sites of the waves' interaction (Fig. [8j;-e). Distribution of 
precipitate (Fig. [8^) inside S approximates skeleton of S (Fig. [Tj:). 
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(e) t = 20 (f) t = 28 



Figure 6. Approximation of Voronoi diagram 03(O) of planar 
shapes O in Voronoi automaton V. Voronoi cell of V in precip- 
itate state are black, cells in refractory states are blank. Excited 
cells are red (gray), (a) representation of O, (b)-(f) development 
of V. 23(0) is represented by black pixels in (f). 
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(a) data 



(b) t = 1 



(c) t = 5 



Figure 7. Schematic representation of data objects and their 
Voronoi diagram or skeleton for precipitation-based approximation 
shown in (a) Fig. [3| (b) Fig. [6j and (c) Fig. § 

6. Summary 

We demonstrated that is possible to solve some problems of computational ge- 
ometry on a discrete model of irregularly arrangements non-uniform vesicles filled 
with chemical mixtures [6l [3 [27] . We shown that a threshold of relative local exci- 
tation density can be used to parameterize excitation and precipitation dynamics, 
enhance results of the computation and reduce noise-to-signal ratio. We proved that 
the same optimal threshold of precipitation works well not only for approximating 
Voronoi diagram of planar sets but also for arbitrary shapes and skeletonisation of 
planar contours. In terms of automata-network based computation we advanced 
our previous results on excitation in Delaunay automata [4] and ^-skeletons [8] . 
The parameterisation developed could be used in the designing experimental labo- 
ratory prototypes of fine-grained compartmentalised excitable chemical processors, 
e.g. using micro-emulsion approach |28j . and developing nano-scale massively par- 
allel computers [TTj . 
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(a) data 



(b) t = 1 




(c) t = 5 



(d) t = 13 




-mm 



(e) t = 26 



Figure 8. Approximation of a skeleton of a star-contour S by V. 
Voronoi cell of V in precipitate state are black, cells in refractory 
states are blank. Excited cells are red (gray), (a) representation 
of a star-contour, (b)-(e) development of V. 
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